H3 in Diffuse Interstellar Clouds: A Tracer for the Cosmic-Ray 

Ionization Rate 

Nick Indriolo 1 , Thomas R. Geballe 2 , Takeshi Oka 3 , and Benjamin J. McCall 1 



ABSTRACT 

Using high resolution infrared spectroscopy we have surveyed twenty sight- 
lines for H^" absorption. Hjj~ is detected in eight diffuse cloud sightlines with 
column densities varying from 0.6 x 10 14 cm -2 to 3.9 x 10 14 cm -2 . This brings to 
fourteen the total number of diffuse cloud sightlines where Hjj" has been detected. 
These detections are mostly along sightlines concentrated in the Galactic plane, 
but well dispersed in Galactic longitude. The results imply that abundant H3 is 
common in the diffuse interstellar medium. Because of the simple chemistry asso- 
ciated with Hjj~ production and destruction, these column density measurements 
can be used in concert with various other data to infer the primary cosmic-ray 
ionization rate, ( p . Values range from 0.5 x 10~ 16 s _1 to 3 x 10~ 16 s _1 with an 
average of 2 x 10~ 16 s _1 . Where H^" is not detected the upper limits on the 
ionization rate are consistent with this range. The average value of ( p is about 
an order of magnitude larger than both the canonical rate and rates previously 
reported by other groups using measurements of OH and HD. The discrepancy is 
most likely due to inaccurate measurements of rate constants and the omission of 
effects which were unknown when those studies were performed. We believe that 
the observed column density of H^" is the most direct tracer for the cosmic-ray 
ionization rate due to its simple chemistry. Recent models of diffuse cloud chem- 
istry require cosmic-ray ionization rates on the order of 10~ 16 s _1 to reproduce 
observed abundances of various atomic and molecular species, in rough accord 
with our observational findings. 



Subject headings: astrochemistry - cosmic rays - ISM: clouds - ISM: molecules 



department of Astronomy and Department of Chemistry, University of Illinois at Urbana-Champaign, 
Urbana, IL 61801; nindrio2@uiuc.edu, bjmccall@uiuc.edu 

2 Gemini Observatory, 670 North A'ohoku Place, Hilo, HI 96720 

3 Department of Astronomy & Astrophysics and Department of Chemistry, University of Chicago, Chicago, 
IL 60637 



-2- 



INTRODUCTION 



I n the past several years Hj~ has been detected in diffuse interstellar clouds (IMcCall et al. 



19981 . 120021 ; iGeballe et al.l Il999l ) where it had been expected to exist in abundances below 
observable limits. This surprising result raised various questions about the diffuse cloud 
environment. The rather simple chemistry of Hjj" allows for only three variable parameters 
in determining its abundance when the steady state approximation is used: the Hj|"-electron 
recombinati on rate, the electron to hydrogen ratio, and the cosmic-ray ionization rate. Pre- 
vious work (iMcCall et al.l 120031 . 120041 ; ICardelli et al.l Il996l ) has shown that the first two of 
these are relatively well constrained. This leaves the cosmic-ray ionization rate as an un- 
constrained parameter. Because the low energy cosmic-rays responsible for most of the 
ionization in diffuse clouds cannot be directly measured in th e solar system , we m ust rely on 



molecules to act as tracers of the ionization rate. Using Hj~, IMcCall et al.l (120031 ) found the 
cosmic-ray ionization rate of molecular hydrogen, £ 2 , to be much larger along the sightline 
to ( Per than the canonical value of ~ 3 x 10~ 17 s -1 . 

Prior to the detection of Hjj~ in diffuse clouds, OH and HD were the molecules of choice 
for estimating the cosmic-ray ionization rate there. Estimates using these molecules required 
determini ng rate constants and modeling various reactions on the pathways to formi ng OH 
and HD (black fc Dalgarnclll977l : iFederman eiTall 1 19961 : lO'Donnell fc Watson! [19741 ) . The 



derived values of the ionization rate tended to agree with the canonical value of ( p , the 
primary cosmic-ray ionization rate, but differ greatly from the value derived from the recent 
H^" measurement toward ( Per (the relation between ( p and (2 is explained in § 4.2 and 
quantified by equation (flOl) ). Of the three molec ules, the simple chemistry of H^" provides the 
most direct determination of ( p (lDalgarndl2006l ). suggesting that measurements of Bj~ should 
produce more accurate results and be a more reliable tracer of the cosmic-ray ionization rate 
than OH or HD. 



The higher ionization rate found by IMcCall et al.l (120031 ) towards ( Per implies the pro 



duction of more H3 , and if generally applicable, could account for th e higher than expected 



column densities found in several diffuse clouds (IMcCall et al.l 120021 ) . However, prior to the 
present work the enhanced ionization rate was known to exist for certain only along one line 
of sight, and thus could have been considered an anomaly. To test if an enhanced ionization 
rate is a general property of the diffuse interstellar medium (ISM), we have performed a 
survey of Hjj" in nineteen diffuse cloud sightlines. Hjj" is detected in eight of the clouds and 
the overall results, including analysis of previous observations by our group, support a higher 
ionization rate. When coupled with further arguments, this strongly suggests that a greatly 
enhanced ionization rate is a typical property of the diffuse ISM. 
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METHODS 



2.1. Observations 



All observations were made using the CGS4 spectrometer ([Mountain et al.lll990l ) on the 
United Kingdom Infrared Telescope (UKIRT) at Mauna Kea. The spectrometer was used 
with its echelle grating, 0.6 arcsec wide slit, and long camera to provide a resolving power of 
40000. Observations were taken in an ABBA pattern in which the target star is nodded along 
the slit so that the spectral image falls alternately on different rows of the array. Suitable 
standard stars were observed throughout each night to account for changing atmospheric 
conditions and air mass. With the exception of the Red Rectangle where reddening is due to 
the ejected envelope of a post asymptotic giant branch star, targets were chosen primarily 
by three criteria: (1) sightlines known to pass through diffuse clouds; (2) early spectral 
type; and (3) bright L-band magnitude. The complete dataset consists of twenty (nineteen 
diffuse cloud) sightlines which were observed intermittently between May 2001 and March 
2005. Relevant information concerning these observations i s given in Table [TJ Two of the 
sightlines from the dataset were examined in iMcCall et al.l (120021 ): HD 20041 and £ Oph. 



However, both have previously been studied using only the R(l, l) 1 transition at 37154.8 
A (vacuum wavelengths are used throughout this paper), whereas the new data coyer the 
R (l, l) u and R( 1 0) tra nsitions at 36680.8 A and 36685.2 A, respectively (see lMcCalll (boooh 



or 



McCall fc Oka (120001 ) for a c omplete description o f the transition notation associated with 



Hg). £ Per was investigated in IMcCall et al.l (120031 ). but here data from three more nights 
of observations are included. Overall, we present twenty sightlines for which new or refined 
Hg column densities or upper limits are calculated. 



2.2. Data Reduction 

The reduction process involves multiple steps and software packages. First, raw data are 
run through Starlink's ORAC-DR0 pipeline which processes UKIRT data. Images are then 
transformed to fits format using Starlink's FIGARol package, and neighboring images are 
subtracted from each other in NOAO's IRAffl package. This subtraction serves to eliminate 
atmospheric background and detector bias levels from the image. Still using IRAF, spectra 



1 http://www. oracdr.org/ 

2 http:/ /www. starlink.rl.ac.uk/star/docs/sun86.htx/sun86. html 
3 http://iraf.noao.edu/ 
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are extracted from the neighbor-subtracted images with the apall routine. These spe ctra are 
imported to IGOR Prc@ where we have macros set up to complete the reduction (jMcCalll 



20011 ). During this import, ripples in the spectrum caused by guiding and seeing fluctuations 
while the CGS4 array is shifted along the spectral direction in successive steps of one-third 
of a pixel during observing are removed. Spectra for each object and standard star are then 
coadded. Some frames may be excluded in this step if they are noisier than normal or happen 
to have cosmic-ray hits on the rows where the spectrum is located. Objects are then ratioed 
with standards to remove atmospheric absorption lines and the continuum level is set to 
unity. The ratioing is an interactive process where the user may vary the intensity scaling 
and/or shift the spectrum in wavelength in order to most effectively remove atmospheric lines 
and obtain the most reliable ratioed spectrum at the wavelengths where the Hjj~ lines are 
expected to appear. A fringing pattern caused by the circular variable order-blocking filter 
is another artifact of CGS4 that needs to be removed. This is accomplished by transforming 
a spectrum into Fourier space using an IGOR macro. The user can then find the peak 
caused by the fringing pattern and interpolate across it. Once the peak is removed, an 
inverse Fourier transform is performed to produce a s pectrum wher e the fringing pattern 



is absent. This method is described in more depth in iMcCalll (120011 ). After removing the 
fringing pattern, spectra are wavelength calibrated using the vacuum rest wavelengths of 
the atmospheric lines. The accuracy of the wavelength calibration is typically ±2 km s _1 . 
Finally, the Hjj" lines are fit with Gaussians and the equivalent widths, column densities, and 
radial velocities are derived. 

For targets with observations on multiple nights, the reduction process above is followed 
through the wavelength calibration step. At that point IRAF's rvcorrect routine is used to 
calculate the velocity of the Earth along a given sightline in the local standard of rest (LSR) 
frame. Each spectrum is then shifted to be in the LSR frame. This puts the interstellar 
absorption lines at the same wavelength for any given date and allows spectra from different 
nights to be coadded. While each spectrum has its continuum set to unity, they originally 
had different exposure times and intensities. To weight the final coadded spectrum properly, 
each spectrum is scaled by its original coadded intensity before being added to the spectra 
from other nights. Once the final multiple night coadded spectrum is produced, it is divided 
by the sum of the scaling factors to reset the continuum level to unity. After this process, 
line parameters are again extracted by fitting the absorption lines. 



4 http: / / www.wavemetrics.com/ 
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2.3. Atmospheric Interference 

Spectra covering the Hjj" doublet are adversely affected by an atmospheric CH4 line 
complex which lies just shortward of 36680.8 A, the rest wavelength of the R(l, l) u line. 
If the spectrum of the target is blue-shifted, the overlap with this line can be significant 
and hinder detection or estimates of line strength. Typical examples of spectra before the 
atmospheric lines are removed via ratioing with a standard star are shown in Figure [TJ 
The CH4 feature here absorbs roughly 50% of the incoming light at 36675.3 A in both the 
standard star (3 Per (middle spectrum) and the object ( Per (bottom spectrum). A weak 
telluric HDO absorption line near 36681 A that can vary with both time and airmass further 
complicates the reduction, especially when the water column density above the telescope is 
high and unstable. When the two spectra are ratioed following the methods in § 2.2, the 
top spectrum in Figure [1] is produced. The two lines show the rest wavelength positions 
of the R(l, l) u line at 36680.8 A and the #(1,0) line at 36685.2 A. Two arrows mark the 
expected positions of the Hjj" lines due the Earth's orbital motion and the radial velocity of 
the absorbing cloud along the line of sight towards ( Per. Clearly the absorption lines 
are much weaker than the atmospheric absorption lines and are barely visible at this scaling 
factor. This illustrates why a careful multi-step reduction process is necessary to detect FPj". 

3. RESULTS 

3.1. Positive Detections 

The fully reduced spectra are shown in Figures [2]l5l Figures [2] and [3] contain spectra from 
sightlines with positive Hjj" detections, and Figures H] and [5] show spectra from sightlines with 
no Hjj" detections or marginal detections. Arrows indicate the position of the Hjj~ doublet 
expected from previous measurements of the gas velocity along each line of sight. These 
velocities are given in Table [2] along with the atomic or molecular species from which they 
were determined. In the figures the typical noise level can be judged by the peak-to-peak 
fluctuations in the continuum well off of the Hjj" lines. Within several Angstroms of 36675.3 A, 
however, the noise in the reduced spectra is several times larger than elsewhere due to the 
presence of the strong complex of CH 4 lines, which both reduces the atmospheric transmission 
and emits excess background. 

The three spectra in Figure [2] have the strongest Hjj~ absorption. In each case the 
lines match up with the arrows so we are confident in the H^" detection. In HD 20041 the 
R(1,0) transition is clearly visible but the R(l,l) u transition may have been affected by 
the aforementioned CH4 line. Both absorption lines are clear and strong in HD 229059. 
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W40 IRS la presents a somewhat confusing case because there may be either one or two 
velocity components. If there is one component, then the two absorption lines to longer 
wavelengths represent the H^ doublet and the shorter wavelength line at 36677 A is a noise 
artifact. If there are two absorption components, then the R(l, 0) transition from the shorter 
wavelength doublet overlaps the R(l, l) n transition from the longer wavelength doublet. 
However, the W40 IRS la spectrum is composed of data from two nights and on only one 
of those nights does the 36677 A feature appear. For this reason we conclude that only the 
lo nger wavelength doublet is real. Our conclusion is consistent with the observations reported 



by lCrutcher fc Chul (119821 ) in which only one velocity component at about 8 k m s 1 LSR was 



seen i n 13 CO emission, but does not match the 2±2 km s 1 LSR reported by IShuping et al. 
( 119991 ) from 12 CO absorption. Note, however, that if the other line is real then it would 



correspond to a cloud with a radial velocity of about -30 km s 1 LSR, which disagrees badly 
with both CO measurements. 

In Figure [3] the absorption lines are again aligned with the arrows marking previously 
observed gas velocities. ( Per and HD 21389 have the highest signal to noise ratios (SNR) 
and their H3" lines are easily identified. The doublet in X Per is relatively clear, but the 
velocities found by fitting the individual line profiles differ by 2.5 km s -1 . Most likely this 
is due to noise affecting the absorption feature. In HD 169454 the R(l, 0) transition of the 
main doublet is clear, but the R(l, l) u line is rather shallow. As in the case of HD 20041 
above, this may be caused by interference from the telluric CH4 line. However, the velocity 
we derive for the R(1,0) transition differs from the previously measured cloud velocity by 
about 3.2 km s _1 . Because of the different velocities and the lack of a clear R(l, l) u line, 
we are not as confident as in the previous cases that these features represent Hg absorption, 
but still consider it a positive detection. The shorter arrows above the HD 169454 spectrum 
mark the e xpected location of the Hj~ d oublet for a high velocity component at 90 km s _1 



reported by lFederman fc Lambertl (119921 ). Our spectrum may indicate absorption features at 
this velocity, but further integration time is needed to determine if the features are real. The 
absorption features in BD -14 5037 both appear to be double peaked, although the signal-to- 
noise ratio of each double peak is low. It also seems that there are small absorption features 
on the shorter wavelength s houlders of the main features. Two velocity components reported 



by iGredel fc Miinchl (119861 ) match well with the centers of the double peaked features and 



the shoulder features, making the detections more believable. 
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3.2. Negative Detections 

There is no definitive evidence for Hjj~ absorption lines in any spectrum in Figure HI 
with the possible exception of o Per. The arrow marking the R(l, 0) transition of H^ in 
o Per matches rather well with a statistically significant absorption feature in the spectrum. 
However, at the expected wavelength of the R(l, l) u line the absorption is too weak to 
attempt a fit given the noise level. Both £ Per and e Per show some amount of absorption 
near the wavelengths where the Hg" lines are expected, but nothing that can be conclusively 
identified as a detection. The Red Rectangle and ( Oph sightlines both have a very high 
SNR. With these clean spectra and no H3 absorption, it is possible to derive strict upper 
limits. HD 147889, 40 Per, and o Sco all have no significant absorption features close to the 
expected wavelengths. Figure [5] contains two non-detections in spectra with typical noise 
levels (HD 168625 and A Cep) and two with high noise levels (HD 21483 and 62 Tau). 



4. ANALYSIS 

By fitting the absorption lines in the spectra with Gaussians we are able to obtain the 
line of sight velocity, full width at half maximum, equivalent width, and the H3" column 
density and its uncertainty. Values for these parameters along all of the observed sightlines 
are given in Table [3j Using these values in concert with the steady state approximation and 
a few reasonable assumptions allows for the calculation of other physical parameters of the 
diffuse clouds along these lines of sight. 



4.1. Reactions 



Below are the three reactions which describe the dominant creation and destruction 
processes for H3 . Reactions (JTl) and (jSJ) show the formation process, while reaction (j3J) 



shows destruction. 



CR + H 2 -> CR + H+ + e - 



Ho + H 9 



H+ + e" 



- H 3 + + H 
H 2 + H or 3H 



(1) 
(2) 
(3) 



First, H 2 is ionized to produce Hj" and an e lectro n. This ionization is assumed to be 
due to a cosmic-ray because iGlassgold fc Langerl (1 19 74 ) showed that low energy cosmic-rays 
will penetrate diffuse clouds while the X-ray flux is attenuated in a thin layer at the cloud 
exterior. They find that for an X-ray with energy 100 eV the optical depth of the cloud 
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reaches unity at an H2 column density of 2 x 10 19 cm -2 . Because most of our sightlines have 
column densities much larger than this we assume that cosmic-rays are the only ionization 
mechanism operating throughout the majority of the cloud. After being ionized the ion 
reacts with H 2 to p roduce Hj~ and H. The second step is many orders of magnitude faster 
than the first step (IMcCall et al.lll998l ). so the formation rate of is proportional to the 
product of the ionization rate and H 2 density. In diffuse clouds the primary channel of 
destruction is electron recombination, which results in either three H atoms or one H atom 
and one H 2 molecule. The destruction rate is given by a rate constant times the product of 
the number densities of H^" and electrons. 



4.2. Calculations 



The steady state approximation assumes that the formation and des truction rates of 
are equal. This approximation yields the equation ( Geballe et al.|[l999"l ) 



n(H 2 )C 2 = A; e n(H 3 +)n(e) (4) 

where n(X) is the number density of species X, ( 2 is the ionization rate of H 2 , and k e is 
the Hjj"-electron recombination rate c onstant. In c ontrast to the steady state approximation, 



time dependent models developed by iLisztl (120071 ) showed that the abundance of is only 



weakly dependent on the cosmic-ray ionization rate when a cloud is young. This age is 
quantified by the ratio of molecular hydrogen to atomic hydrogen n(H 2 )/n(H) where smaller 
values correspond to younger clouds. The abundance becomes weakly dependent on the 
cosmic-ray ionization rate when n(H 2 )/n(H) < 0.05. This value corresponds to a molecular 
hydrogen fraction (defined below in equation (Q) of / < 0.09. Because the H 2 fractions in 
all of the clouds we observed are more than double this value, we neglect time dependence 
and use the steady state approximation. 

Assuming that gas is uniformly distributed in the cloud, we can substitute the column 
density divided by the path length for the number density. By doing this and solving for the 
ionization rate we obtain 

L n(H 2 ) 

We further assume that nearly all hydrogen is either in the atomic or molecular state, and 
define the molecular hydrogen fraction (the fraction of hydrogen nuclei in molecular form) 

as 

= 2n(H 2 ) 

1 ~ 71(H) + 2n(H 2 ) 1 ' 
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where the denominator is the number density of hydrogen nuclei, n^. Solving for n(R.2) and 
plugging the result back into the ionization rate equation we find 



C 2 = N(E+) 



k P 2 n(e) 



(7) 



3 ' L f n R 

In this form it is possible to measure or estimate all of the variables on the right hand side 
of the equation, so we can derive values for the ionization rate of molecular hydrogen. 

The electron recombination rate constant is given by the equation 

K = -1.3 x 1(T 8 + 1.27 x io- 6 T-°- 48 (cm 3 s" 1 ) (8) 



from iMcCall et al.l (120041 ) which is valid when the electron temperature, T e , is between 10 
K and 4000 K. While T e is not directly measured, it can be approximated by the excitation 
temperature derived from the J = and J = 1 levels of molecular hydrogen, Tq\. This 
temperature is calculated from measurements of the column densities of the two levels. In 
sightlines without these measurements we adopt a value of 60 K. The fact that the observed 
lines of the J = and J = 1 levels of H2 are saturated indicates that few photons are present 
in the interior of diffuse clouds to radiatively pump these levels. This means that collisions 
will dominate the equilibrium between these levels. However, the J = and J = 1 levels of 
H2 have different nuclear spin configurations an d thus require collisions with species such as 
H + or Hjj~ to interconvert (ISnow fc McCallll2006l ). If collisions with protons are the dominant 
factor in determining the relative population of the J = and J = 1 levels though, then 
Tqi should represe nt the proton kinet ic temperature and thus the kinetic temperature of 
the gas in general (ISavage et al.lll977l ). While electrons produced by photoionization may 
begin with much hig her temperatures, they should thermalize quickly via collisions with H2 
( IMcCall et al.l 120021 ) . Because T e and T 01 should both nearly equal the kinetic temperature 
of the gas, we substitute T 01 for T e in equation (jSJ). 

Assuming that nearly all electrons in diffuse clouds are produce d via the ionization of 
C to C + and that nearly all atomic carbon has been singly ionized (Ivan Dishoeck & Black 
1986 ), the carbon to hydrogen ratio should approximate the electron to hydrogen ratio. 
Cardelli et al.l (119961 ) found this value to be about 1.4 x 10~ 4 in multiple diffuse clouds. 



Because of the relative uniformity of this ratio in all six of their sightlines, we adopt a single 
average value for use in all of our calculations. 

The molecular hydrogen fraction is dependent on H and H 2 number densities, quantities 
which surely vary through the cloud, but whose variations are not readily measurable. Be- 
cause fluctuations in number density cannot be directly measured, A ye use column densities 
in place of the number densities in equation ([6]) and calculate what ISnow k, McCalll (120061 ) 
refer to as f N in clouds where we have measurements of the H and H 2 column densities. 



-10- 



However, f N is most likely an underestimate of / in the more molecular regions of the cloud 
which contain higher concentrations of hLt. This is because atomic hydrogen is more widely 
distributed t han molecular hydroge n and column densities measure material along the entire 
line of sight (jSnow &: McCalll 120061 ) . Since the measurement of iV(H) includes material not 
associated with H2, / underestimates the H2 fraction in the molecular region. For sightlines 
where measurements are lacking we use / = 0.67, the value for which the column densities 
of H and H 2 are equal. 

When possible, estimates of the number density of hydrogen nuclei, %, are adopted 
from the literature based on v arious atomic and molecular diagnostics. In eleven of our 
sightlines ISonnentrucker et al.l (120071 ) used the observed rotational excitation of C2 to infer 
the sum of the H and H 2 number densities. This was done by comparing models with various 
temperatures and number densities to the measured column densities of all the excited states 
and choosing the best fit. For sightlines where the average value of / is known, they converted 
n(H + H 2 ) to %. For sightlines where / is not kno wn, ou r adop ted value of / = 0.67 is used 
to perform the conversion. In two additional cases lJural (119751 ) measured column densities 
of H and the J = 4 excited level of H2, and with some assumptions estimated the product 
Rur, where R is the rate at which H 2 forms on grains. Adopting a typical value for R 
then allowed for the computation of tir. In one more sightline (40 Per) iJenkins et al.l (119831 ) 
estimated the thermal pressure from measurements of the J = 0, 1, and 2 fine-structure 
levels of C for a kinetic temperature of 80 K. Using this pressure estimate and the H 2 
temperature (T i = 63 K), we calculate n H - Unfortunately, the results obtained for a given 
sightline by using each of these methods can be sig n ifican t ly dif ferent. For example, in the 
sightline toward ( Oph ISonnentrucker et al.l (120071 ). lJural (119751 ). and IJenkins et al.l (119831 ) 
derived values of 215 cm -3 , 90 cm -3 , and 117 cm -3 , respectively for %• Because of the 
uncertainties involved with each method and the different final results, the number densities 
we use are probably uncertain by about a factor of two. For cases where no number density 
has been determined, a value of 250 cm -3 is adopted. 



Again assuming a uniform distribution of gas in each cloud, we divide the total hydrogen 
column densities by these number densities to obtain pathlengths: 

Nh = N(H) + 2iV(H 2 ) 
n H n H 

In sightlines where no H and H 2 column densities have been determined, the total hydrogen 
column de nsity is estimated from the color excess a nd the relation Nh rs E(B — V) x 5.8 x 



(9) 



10 cm (jBohlin et al.lll978l : lRachford et al.ll2002l ). As the pathlength is calculated directly 



from the hydrogen number density, it is also uncertain by about a factor of two. Using the 
above relations and approximations, we calculate the ionization rate of molecular hydrogen, 
( 2 - However, most studies examine the primary cosmic-ray ionization rate per hydrogen 
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atom, (p. While the ionization efficiencies of H and H2 are dep endent on factors suc h as 



the helium abundance and ratio of molecular to atomic hydrogen (IDalgarno et al.lll9991). we 



adopt a more simplified approach and use the conversion factor given by iGlassgold fc Langer 



(|l974h : 

C2 = 2.3C P . (10) 

This conversion stems from the fact that H 2 contains two hydrogen atoms, so the ionization 
rate is nearly twice as high. Also, ( p only accounts for the initial (primary) cosmic-ray 
ionization while (2 includes ionization from energetic secondary electrons which were created 
in the first ionization event. With equation (flOl) we convert our values of (2 to ( p so that 
they can be directly compared to previous observations. The resulting values for the primary 
cosmic-ray ionization rate and the specific estimates used for each sightline are shown in 
Table HI For comp leteness, we have performed the same analysis for ten sightlines from 



McCall et al.l (120021 ) and also included the results in Table HI 



With all of the assumptions we have made, it is important to investigate the uncertainties 
that will propagate to the cosmic-ray ionization rate. Substituting equations Qj and ffTUl) 
into equation ([7]) gives the primary cosmic-ray ionization rate as 



n(e) 
n H 



In this equation, the molecular hydrogen fraction, /, and the number density of hydrogen 
nuclei, are the two most uncertain parameters (note that the % in the denominator is 
part of the ratio n(e)/% which is well determined). Because ( p is directly proportional to 
%, any increase or decrease in produces a corresponding increase or decrease in ( p . As 
previously mentioned, the uncertainty in nn is probably about a factor of 2. On the other 
hand, ( p is inversely related to /. The H2 fraction is by definition between zero and one, 
and most of our measured values are around 0.5. For sightlines with the adopted value of 
0.67, the maximum increase is a factor of 1.5. We take this factor to be an approximation 
for the uncertainty in /. Because f N should always be an underestimate of / in measured 
sightlines, we only consider increasing the molecular hydrogen fraction for those sightlines. 
Taking into account the uncertainties in both «h and /, the true value of ( p in sightlines with 
measurements of f N is likely between one third and twice our derived estimate of ( p . For 
sightlines with no measured H2 fraction, we allow / to vary both up or down by a factor of 
1.5. This results in a possible cosmic-ray ionization rate between ( p /3 and 3( p . These limits 
arise when the most extreme variations in both / and % are substituted into equation (TIT]) . 
However, / tends to be higher when % is higher because the rate of H 2 formation scales as 
the square of %. This suggests that it is probable that / and n-^ will vary in the same way, 



so that the above analysis most likely overestimates the range of possible ionization rates. 



4.3. H3 Temperature 



Another property of the gas to be examined is the excitation temperature, determined 
from relative populations of the different rotational states of H3 . For the two lines we have 
observed the temperature may be determined from the equation 



N, 



N, 



oil ho 9 ortho ^— l\E / kT 



para 



2e -32.87/T 



Qpara 



(12) 



taken from iMcCall et al.l (119981 ). In this case ortho refers to the population of the (1,0) 
state and para to the (1,1) state. The g's are statistical weights, AE is the energy difference 
between the states, k is the Boltzmann constant, and T is the temperature. Excitation 
temperatures derived from this equation are shown in Table HI 

If the rotational (de-)excitation of is dominated by collisions with H2, then the tem- 
perature measurements from both species should be similar. However, this is not the case. 
Most H 2 temperatures are around 60 K while the Hj~ tempera tures are typically about 30 K. 
This same discrepa ncy was described by IMcCall et al.l (120031 ). In their model calculation of 



thermalization, lOka &: Eppl (120041 ) have shown that the (1,1)/(1,0) excitation tempera- 



ture is always lower than the H2 temperature because of cooling by fast spontaneous emission 
from the (2,2) to (1,1) s tate. For the typic al cloud conditions in this paper (% ~ 250 cm -3 , 
T ~ 60 K) the model of lOka fc Eppl (120041 ) produces an H3 excitation temperature of about 
50-55 K which is significantly higher than the observed values of about 30 K. The source of 
this discrepancy remains unclear. 



5. Discussion 

5.1. Inferred Ionization Rates 

Values of the cosmic-ray ionization rate for diffuse clouds observed here as well as those 
determined for other diffuse clouds observed previously by us are given in the right hand 
column of Table HI The detected values in the lines of sight to fourteen sources cover the 
range 0.5-3.2 xlO 16 s _1 . Upper limits, which are given for fifteen diffuse clouds, are consistent 
with this range of ionization rates, with the possible exception of HD 168607. While most of 
the detections of are confined to the Galactic plane, they are widely dispersed in Galactic 
longitude. We therefore conclude that the values of the cosmic-ray ionization rate listed in 



13 



Table 0] are typical for Galactic diffuse interstellar clouds. 



A few of the sightlines we investigated have been studied previously to derive cosmic-ray 
ionization rates. All of these studies used column densities of either OH, HD, or both in 
their calculations. Because the formation pathways of OH and HD include the ionization 
of atom ic hydrogen, they can b e used to determine the H ionization rate. Most of these 



studi es (IBlack fc Dalgarnol 1 19771 ; iBlack et al.l 1 19781 ; lHartquist et al.l Il978al ; iFederman et al. 



1996 ) then derived the prim ary cosmic-ray ionization rate from the H ionization rate, but 
O'Donnell fc Watsonl (11974 ) did not because they still considered ionization via X-rays to 
be important. Our values of the primary ionization rate for ( Per, o Per, e Per, £ Per, and 
( Oph are shown in Table [5] along with the rates derived from OH and HD measurements as 
well as clou d modeling. For ( Per o ur va lue is over an order of m agnitude larger than those 
reported by lHartquist et al.l (ll978bl ) and IFederman et al.l (119961 ) . While the rest of our new 
measurements in Table [5] are only upper limits, these are also typically orders of magnitude 
larger than previously published values. The only exception is o Per where both papers cite 
values of ( p about one fourth to one half our upper limit. 



Various model calculations were performed by Ivan Dishoeck &: Black! (119861 ) to inves- 
tigate three of the sightlines that we study here: ( Oph, ( Per, and o Per. In creating 
these models they used the most recent measurements of rate constants and the column 
densities of diagnostic species such as H and H2 as input parameters. By varying a few 
uncertain parameters, they would then generate lists of predicted column densities for many 
atomic and molecular species under slightly different conditions. When their paper was 
written, it was believed that the Hj~-electron reco mbination rate constant was much lower 



than the currently accepted value. iSmith fc Adamsl (ll984j ) repor ted an upper limit corresp- 
nding to 10~ 7 cm 3 s _1 at T = 40 K, and lAdams fc Smith! (119871 ) lowered the upper limit to 
10~ n cm 3 s" 1 at T = 80 K. D ue to the wide range of possible recombination rate constants, 
van Dishoeck &: Black! (119861 ) performed calculations using both 10~ 7 and 10~ 



-10 cm 3 g 1^ 



The cosmic-ray ionization rates from their paper listed in Table [5] were computed by deter- 
mining ( p necessary to reproduce observed OH column densities when k e = 10 -7 cm 3 s -1 . 
We choose to compare these ionization rates to ours because we obtain k e = 1.6 x 10~ 7 
cm 3 s -1 w hen T = 60 K is used a s the i nput temperature in equation (jSj). The value of ( p 
inferred by Ivan Dishoeck fc Black! (119861 ) is about the same as ours for £ Per, but the lower 
limits they derived for ( Oph and o Per are larger than our upper limits for both of those 
sightlines. 



For their models that used k e = 10~ 10 cm 3 s _1 , Ivan Dishoeck &: Black! (119861 ) obtained 
cosmic-ray ionization rates that are about a factor of 1 to 5 times smaller than ours. From 
these models they also predicted the column density of H^ along each sightline. Their results 
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are all on the order of iV(H3 ) ~ 10 cm , which is a few times la rger than the observed 



colum n densities or upper limits in any of these sightlines. Because Ivan Dishoeck fe Black 
(119861 ) use only a slightly smaller cosmic-ray ionization rate (corresponding to the formation 
rate) but a much smaller recombination rate (corresponding to the destruction rate), their 
prediction of an Hg column density similar to observed values seems somewhat serendipitous. 
In addition to the overestimate of the Hjj" colu mn density, a sma ll Hjj~-electron recombination 
rate constant may have further consequences. iDalgarnd (120061 ) noted that a small value of 
k e may have been responsible for underestimates of the primary cosmic-ray ionization rate 
in the past. This is because a slower destruction rate requires a slower formation rate to 
produce a given abundance. 

In addition to the slow recombination rate, there are some other possible explanations for 
differences be tween the cosm i c-ray ionization rate inferred from H^" and those inferred from 
OH and HD. iLe Petit et al.l (120041 ) pointed out that the rate constant associated with the 
endothermic charge transfer from H + to O varies over the temperatures typically associated 
with diffuse clouds. This means that the OH producti on rate is temperature de pendent. The 
ionization rates towards ( Per and ( Oph quoted in lHartquist et al.l (jl978bl ) were derived 
using temperatures of 120 K and 110 K, respective l y, for the warm compon ents of the 
cloud models along each sightline (IBlack et al.l Il978l ; iBlack fc Dalgarnol 119771 ). As these 
temperatures are about twice as large as the values determined from H 2 , their OH production 
is much more efficient. The result is a smaller cosmic-ray ionization rate needed to produce 
the observed OH column density tha n if a lower temperatur e had been used. This problem 
was addressed by the later models of Ivan Dishoeck fc Blackl (119861 ) where tempertaure and 
density were varied as functions of cloud depth. 



Le Petit et al.l (120041 ) went on to make a comprehensive chemical model of the cloud 



towards ( Per. They determined the value of ( p that would best reproduce all observed 
atomic and molecular column densities to be 2.5 x 10~ 16 s -1 , which is in good agreement 
with our estimate of 3.2 x 10~ 16 s _1 . Th e difference in these y alues may arise because we 
assume a uniform distribution of gas while ILe Petit et al.l (120041 ) invoke a three phase model 
which includes diffuse gas, dense gas, and magnetohydrodynamic shocks. 

Two different effects may lead to underestimates of ( p from measurements of HD. The 
first has to do with an overestimate of the total deuterium to hydrogen ratio no/%. This 
ratio can be used to estimate the molecular deuterium fraction, iV(HD)/iV(H2). However, 
the observed values of A^(HD)/A^(H 2 ) are about an order of m agnitude smaller than those 
predicted by n^/riu- To explain this discrepancy, iLisztl (120061 ) argued that the atomic deu- 
terium fraction must be larger than the total deuterium fraction. This means that ap- 
proximating hd/wh with iV(D)/iV(H) overestimates the total deuterium to hydrogen ratio. 
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Federman et al.l (119961 ) showed that the cosmic-ray ionization rate is inversely related to the 
deuterium fraction, so an overestimate of n-£,/nn will underestimate ( p . 



Second l y. iLisztl (120031 ) emphasized the importance of grain neutralization proposed by 
Lepp et al. ( 19881 ). This process reduces the number of H + ions in the gas through charge 
transfer with small grains. By lowering the abundance of H + , the production rate of HD will 
decrease. This is because HD formation is dependent upon the reaction involving the charge 
transfer from H + to D. Since neutralization slows down HD production, a larger value of ( p 



i s nee ded to create a given abundance than if the effect were not taken into account. [Liszt 



( 120031 ) used a model which includes grain neutralization and showed that both Hjj~ and HD 



column densities can be reproduced with a single ionization rate of ( p > 2 x 10 



- 16 s- 1 . 



Since 



OH formation is dependent on a similar charge transfer reaction, grain neutralization and 
thus a larger cosmic-ray ionization rate may be necessary in its analysis as well. 



McCall et al.l (120031 ) studied Hjj~ in the sightline towards ( Per. Using nearly the same 



analysis as this paper, they inferred a value of C2 = 1-2 X 10 s which is equivalent to 
( p = 5.2 x 10~ 16 s _1 shown in TableEJ This higher ionization rate i s due to various diffe rences 
in input parameters. In terms of the parameters in this paper, iMcCall et al.l (120031 ) used 
1.5/c e , 1.2A^(H^), 1.2riH, and 0.8n(e)/n-n for the following reasons. The H^"-electron recom- 
bination rate constant differs because they approximated the electron temperature with the 
temperature instead of the H2 temperature in equation (jSj). Further observations have 
more than doubled the total integration time so that t he spectrum and Ht column density 
change slightly between papers. The value of n H used by lMcCall et al.l (120031 ) was an average 
number density computed fro m various measurements, w hereas the value used in this paper 
comes from the C2 analysis of ISonnentrucker et al.l (120071 ). Finally, we have adopted a single 
value of n(e)/n H to be used in all calculations while they used H 2 and C + column densities 
measured towards ( Per. 

While all of the observations and models above are viable methods for finding the 
cosmic-ray ionization rate, we believe that the use of H3 should produce the best results 
due to its relatively simple chemistry. Using either OH or HD to calculate ( p requires more 
measurements, more assumptions, and more variable parameters than using H^. More pa- 
rameters give the opportunity for a greater uncertainty to accumulate during the calculation. 
Fewer uncertainties coupled with advances in instrumentation lead us to speculate that the 
cosmic-ray ionization rates inferred from H^ may be the most accurate to date for diffuse 
clouds. However, improved estimates of / and n H , the two most uncertain values in our 
calculations, would make Hjj" an even better probe of the cosmic-ray ionization rate. 
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5.2. Theoretical Ionization Rates 



Several theoretical calculat i ons of ( p have been perfo r med in the las t half-century (IHayakawa et al 



196ll ; ISpitzer &: Tomaskolll968l ; iNath &: Biermannl 11994 ; IWebberlll998l ). In these papers the 
authors derived a cosmic-ray ionization rate starting from the observed flux of cosmic-rays in 
our solar system. Unfortunately, there are large uncertainties associated with this method. 
The cosmic-ray spectrum is well measured above about 1 GeV, but lower energy particles 
are deflected from the inner solar system by the magnetic field coupled to the solar wind. 
The particles which are most important for ionizing species in diffuse clouds are likely those 
with energies from about 2 to 10 MeV. Since this portion of the spectrum cannot be di- 
rectly measured, the flu x at low energies must somehow be extrapolated from existing data. 
Hayakawa et al.l (ll96ll ) assumed that the power law which applies to the flux of high en- 
ergy cosmic-rays continues down to 10 MeV where the spectrum peaks and then decreases 
linearly with energy . From these assumptions they derived an ionization rate of 10 -15 s _1 . 
Spitzer fc Tomaskd (119681 ). however, fit a curve to measurements of cosmic-rays with energies 
near 100 MeV that also matches the high energy spectrum power law. With this method, 
their spectrum peaks around 100 MeV and falls off for lower energies. The result of using 
their fit is a lower limit of 6.8 x 10~ 18 s _1 . In the same paper they derived an upper limit of 
1.2 x 10 -15 s _1 v ia arguments th at low energy cosmic-ray protons are accelerated in Type I 
supernova shells. IWebberl (119981 ) used data from the Pioneer and Voyager spacecraft as they 
travelled outward in the solar system where the weaker solar wind allows for the detection 
of lower-energy cosmic-rays. These data were then combined with previous observations 
to infer the interstellar proton spectrum. Using this proton spectrum and a heavy nuclei 
spectru m both with lo w energy cut-offs at 10 MeV and an electron spectrum cut-off below 
2 MeV. IWebberl (119981 ) calculated the primary cosmic-ray ionization rate to be (3 — 4) x 10~ 17 



s _1 . Our ionization rates fall neatly within the bounds formed by these studies and so are 
not inconsistent with constraints based on direct cosmic-ray measurements and theoretical 
particle physics. 



5.3. The Ionization Rate in Dense Clouds 



In contrast to our findings in diffuse clouds, the cosmic-ray ionization rate in dense 
clouds does seem to agree with the canonical value. Observations of Hg towar ds dense clouds 
have found column densi ties roughly the same as those seen in diffuse clouds ( jGeballe fc Oka 
19961 ; iMcCall et al.lll999l ). These measurements have been used to calculate the product ( 2 L. 
When ( 2 is taken to be the canonical value of ~ 3 x 10~ 17 s _1 , the resulting pathlength is on 
the order of a parsec. This is a typical size for dense clouds as measured by other methods 
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such as extinction mapping. Since should be a reliable tracer for the cosmic-ray ionization 
rate in both environments, there must be some mechanism c ausing the difference bet ween 
dense and diffus e clou ds. One possibility examined by both ISkiUing fe Strong (119761 ) and 



Padoan fc Scald (120051 ) is cosmic-ray self-confinement. In this process cosmic-rays generate 
Alfven waves which can effectively confine the lower energy particles (< 100 MeV) to diffuse 
material, thus preventing them from entering dense clouds. Because cosmic-rays in the 
1-100 MeV range are the most efficient at ionization, self-confinement naturally leads to a 
higher ionization rate in diffuse clouds than in dense clouds. Another possibility is that there 
is a previously unrecognized high flux of low energy cosmic-rays that can penetrate diffuse 
but not dense clouds. Assuming that typical column densities of diffuse clouds are of order 
10 21 cm" 2 and those of dense cloud s are of order 10 23 cm -2 , cosmic-rays with energies ~2- 
20 MeV (ICravens &: Dalgarnolll978l ) would contribute to the ionization rate only in diffuse 
clouds. As we foresee no observational techniques that would distinguish between these 
two possibilities, a resolution to this question will depend on more sophisticated theoretical 
treatments. 



6. SUMMARY & CONCLUSIONS 

We have surveyed twenty sightlines and detected along eight of them. Column den- 
sities are measured for these eight sightlines, and upper limits set for the remaining twelve. 
Besides a concentration near the Galactic plane, there seems to be no clear correlation be- 
tween location in the sky and detecting , so it is unlikely that we are observing anomalous 
regions in the Galaxy. Instead, finding Hg in so many sightlines suggests that it is ubiquitous 
in the diffuse ISM. 

From the Hjj~ column densities and the steady state approximation, we derive cosmic- 
ray ionization rates for the ninete en diffuse cloud sightlines in this study along with ten 



sightlines from iMcCall et all (120021 ). Typical values are on the order of ( p ~ 2 x 10 s 
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which falls within theoretical constraints. While this is an order of magnitude larger than 
most previously inferred values, there are several possible explanations for the discrepancy. 
The most likely candidates are rate constants with uncertain measurements and physical 
or chemical effects not included in past models. Newer models that do take into account 
these factors require cosmic-ray ionization rates very similar to our inferred values. Coupled 
with these models, our widespread detection of Hg in diffuse clouds supports the idea that 
the typical cosmic-ray ionization rate in such regions should be revised upward by about an 
order of magnitude. 
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Fig. 1. — Typical examples of spectra near the Hjj" ortho-para doublet. The spectra in this 
and all other figures have been offset in intensity for clarity. The bottom spectrum is ( Per 
from 2001 September 5. The middle spectrum is the standard star (3 Per from the same date. 
The top spectrum is ( Per ratioed with the standard star. Two arrows show the expected 
location of Hjj" absorption which can barely be seen here due to the scaling. The vertical 
lines are at the rest wavelengths of the H3 lines. 
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Fig. 2. — Spectra showing strong detections of the H;j~ doublet near 36680 A. All spectra 
have been Doppler shifted into the rest frame of the LSR. Arrows show where the lines are 
expected due to previous gas velocity measurements, which are given in Table El 
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Fig. 3. — Same as Figure 2 except showing more typical strength detections of Hjj" 
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Fig. 4. — Same as Figure 2 except showing non-detections. 
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Fig. 5. — Same as Figure 2 except showing non-detections. 
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Table 1. Observations 







Date(s) of 




Integration Time 


Obiect 

J 


HD number 


Observation 


Standard 


(sec) 


COph 


149757 


2001 May 24 


(3 Lib 


576 






2001 May 25 


(3 Lib 


1344 


HD 147889 


147889 


2001 May 24 


5 Sco 


5760 


A Cep 


210839 


2001 May 24 


a Lyr 


4320 


HD 169454 


169454 


2001 May 25 


(3 Lib 


4224 


W40 IRS la 




2001 May 25 


7] Oph 


1800 






2001 May 26 


5 Sco 


1800 


o Sco 


147084 


2001 May 26 


a Lyr 


960 


HD 168625 


168625 


2001 May 26 


i] Oph 


1440 


HD 229059 


229059 


2001 May 27 


a Cyg 


1800 






2001 Sept 5 


a Cyg 


1800 


BD -14 5037 




2001 May 28 


i] Oph 


2016 


HD 20041 


20041 


2001 Sept 5 


(3 Per 


1728 


HD 21389 


21389 


2001 Sept 5 


(3 Per 


1536 






2002 Dec 30 


S Per 


1152 






2002 Dec 31 


(3 Per 


1536 






2003 Jan 1 


7] Tau 


1152 


(Per 


24398 


2001 Sept 5 


f3 Per 


2304 






2002 Dec 30 


f3 Per 


1152 






2002 Dec 31 


(3 Per 


1152 






2003 Jan 1 


S Per 


1440 


o Per 


23180 


2002 Dec 30 


(3 Per 


2304 






2002 Dec 31 


(3 Per 


1152 


f Per 


24912 


2002 Dec 30 


f3 Per 


1920 






2002 Dec 31 


(3 Per 


1920 






2003 Jan 1 


8 Per 


1536 


X Per 


24534 


2002 Dec 31 


(3 Per 


1920 






2004 Jan 22 


HD 17573 


3600 






2005 Jan 5 


HD 17573 


4320 






2005 Jan 6 


i] Aur 


2880 
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Table 1 — Continued 







Date(s) of 




Integration Time 


Object 


HD number 


Observation 


Standard 


(sec) 






2005 Mar 3 


i] Aur 


5040 






2005 Mar 4 


7) Aur 


5040 


62 Tau 


27778 


2003 Jan 1 


(3 Per 


2688 


e Per 


24760 


2004 Jan 23 


r) Aur 


1800 






2005 Jan 5 


i] Aur 


1440 






2005 Jan 6 


i] Aur 


1440 


40 Per 


22951 


2004 Jan 23 


HD 17573 


3600 






2005 Jan 6 


i] Aur 


2880 


Red Rectangle 


44179 


2005 Jan 6 


k Ori 


720 


HD 21483 


21483 


2005 Jan 25 


HD 17573 


5760 
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Table 2. Previous Measurements of ISM Gas Velocities 



Object 


^LSR 

(km s _1 ) 


Species 


Reference 


tt"pi r\r\r\ a -i 

HD 20041 


-1.6 


K I 


1 


HD 21389 


-0.5 


Si T T 

CH 


1 


C Per 


6.9 


f i T T 

CH 


1 


X Per 


6.6 


CH 


1 


TTl — \ 1 f~* r\ A V A 

HD 169454 


5.3 


/- i T T 

CH 


1 




90 


Na I & Ca II 


2 


T T T ~\ c\c\f\r\ V f\ 

HD 229059 


4.04 


K I 


3 


BD -14 5037 


8.2 


c 2 


4 




18.2 


c 2 


4 


W40 IRS la 


8 


1 1nnt 


5 




10.3 


CH 


i 

1 


40 Per 


6.8 


K I 


3 


o Per 


7.3 


CH 


1 


e Per 


2.2 


K I 


6 


f Per 


0.55 


K I 


6 


62 Tau 


4.8 


CH 


1 


Red Rectangle 


-0.8 


K it & I f 


7 


o Sco 


2.29 


K I 


6 


HD 147889 


2.7 


K I 


1 


COph 


-1.0 


CH 


1 


HD 168625 


6 


Ca II 


8 


A Cep 


-1.7 


CH 


1 



References. - (1 ) Welt y (pr ivate communica t ion); (2) iFederman fc Lambertl f 19921 ); 
(3) Ichaffee fc Whit3 Jl982h : (4) lOredel fc Miinri] Jl986h : (5) ICrutcher fc Chul Jl982h : (6) 
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Weltv & Hobbs 


(2001): (7) Hobbs et al. 


(2004 


); (8) 


Rickard 


(1974) 



Note. - wlsr is the velocity of the interstellar gas in the local standard of rest frame. 
Unless noted, all lines were measured in absorption. 

^measured in emission 
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Table 3. Absorption Line Parameters 







^LSR 


FWHM 


W x 




^(H 3 + ) 


/j(7V) 


Object 


Transition 


(km s _1 ) 


(km s _1 ) 


/ 2 \ 

(A) 


(A) 


(10 14 cm' 2 ) 


(10 14 cm 


LlL) ZUU41 


JL^l, 1 J 


1 /i 

-1.4 


8 A 
o.4 


O 01 7 
U.U1 ( 


n nn/i 

U.UU4 


n 7n 
u. i u 


U. 10 




d/i r\\ 
-K(l, U J 


1 t; 
-1.0 


O fi 

y.o 


O OQfi 
U.U00 


n nn/i 
U.UU4 


n oi 
u.y i 


n i n 

U. 1U 


PTT"1 91 'iSQ 

nu ziooy 


D/1 1 Am 
iXl 1, 1 J 


1 Q 

-i.y 


1 1 Q 
1 1 .y 


n nno 
u.uuy 


n nno 
u.uuz 


n Q.Q 
u.oy 


n n7 
u.u / 




D/i ro 

-fl^l, U J 


n 9 
-u.z 


i ^ fi 

10.0 


n m p 

U.U10 


n nno 
u.uuz 


n 4i 

U.41 


n 

u.uo 


C, rer 


d/i 1 A M 
xti i, i ) 


/ .0 


Q 7 

y. / 


ooio 

U.U1U 


n nm 

U.UUl 


n a q 

U.40 


n r\K 
u.uo 




d/i nA 
ix(l, U J 


P Q 

O.o 


8 8 
o.o 


010 
U.U1U 


n nm 

U.UUl 


n op 

U.Z0 


n no 
u.uo 


Y Per 


D/1 1 \u 
ill 1, 1 J 


o.o 


1 1 fi 
11.0 


011 
U.Ull 


n nno 
u.uuz 


n ap 

U.40 


n 1 n 

U. 1U 




d/i ro 

JT^l, U I 


0. ( 


Q P. 

y.o 


019 
U.U1Z 


n nn9 
u.uuz 


n 9.1 

U.Ol 


n nfi 
u.uo 


PPP~> 1 fiQ/1 ^/l 

nu ioy404 


P/1 1 A M 
ill 1, 1 J 


9 fi 
z.0 


111 
11.1 


o ooc; 
u.uuo 


n nn9 
u.uuz 


n 91 

U.Z1 


n na 

U.Uo 




p/i ro 

JT^l, U J 


9 1 
Z. 1 


1 n 8 

lU.o 


014 
U.U14 


n nno 
u.uuz 


n o.^ 
u.oo 


n 

u.uo 


nu zzyuoy 


D/1 1 \M 
iXl 1, 1 J 


^ Q 

o.y 


19 9 
1Z.Z 


0^8 
U.UOo 


n nno 
u.uuo 


9 49 
Z.4Z 


n 1 4 

U. 14 




p/i ro 


A A 
4.4 


1 9 A 
LZA 


0^0 


n nno 
u.uuo 


1 48 
1 .4o 


n no 
u.uy 


RD 1 /I CAQ7 

jdJJ -14 0U0 / 


D/1 1 A M 
xtl 1, 1 ) 


1 7 fi 
1 / .0 


O fi 

y.o 


OOIO 
U.U1U 


n nnQ 
U.UUO 


n /in 
U.4U 


n 1 q 

U. 10 




p/i nA 


1 8 fi 
lo.O 


10/1 
1U.4 


OOQ 

u.uuy 


n nno 
u.uuo 


n 9Q 
u.zo 


n ns 

U.Uo 


W4U ltxD ra 


P/1 1 AM 

JXl i, i J 


/ .0 


10 8 
lU.o 


o oc;i 

U.UOl 


n nns 

U.UUo 


9 19 

z. 1Z 


n qq 
u.oo 




pfi ro 

JT^l, U J 


8 9 
o.Z 


fi fi 
0.0 


u.uou 


n nnfi 
u.uuo 


1 9fi 
1 .ZO 


n 1 fi 

U. 10 


11F> 91 /ICQ 


R/1 1 \ u 
JXyL, 1 J 




1 n 

1U 


^0 OQR 
<U.U00 




^1 i^Q. 
< 1 .00 






d/i n\ 

-fl^l, U ) 




i n 

1U 


"^U.UOO 




^n qq. 
^u.yo 




/in Por 

4U Jrei 


Pf1 A\ u 
JXl 1, 1 J 




1 n 

1U 


^0 01 ^ 
^U.UIO 




^n fin 

*^U.OU 






p/1 f\\ 

-TL^l, U J 




1 n 

1U 


^0 01 
^U.UIO 




<--n Qfi 

^U.OO 




o Per 


D/1 1 A" 
iXl 1, 1 J 




1 n 

1U 


^n nno 
^u.uuy 




^n 49 

^U.4Z 






#(1,0) 




10 


<0.009 




<0.27 




e Per 


#(1,1)" 




10 


<0.009 




<0.39 






#(1,0) 




10 


<0.009 




<0.24 




f Per 


#(l,l) u 




10 


<0.009 




<0.33 






#(1,0) 




10 


<0.009 




<0.21 




62 Tau 


#(1,1) M 




10 


<0.045 




<1.89 






#(1,0) 




10 


<0.045 




<1.14 




Red Rectangle 


#(1,1) M 




10 


<0.006 




<0.21 






#(1,0) 




10 


<0.006 




<0.15 
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Table 3 — Continued 





^LSR 


FWHM 


w x 




iV(H 3 +) a(N) 


Object 


Transition (km s -1 ) 


(km s- 1 ) 


(A) 


(A) 


(10 14 cm- 2 ) (10 14 cm" 2 ) 


o Sco 


R(i,ir 


10 


<0.009 




<0.36 




#(1,0) 


10 


<0.009 




<0.21 


HD 147889 


R(i,ir 


10 


<0.009 




<0.39 




R(1,0) 


10 


<0.009 




<0.24 


COph 


R(l,l) u 


10 


<0.003 




<0.18 




#(1,0) 


10 


<0.003 




<0.12 


HD 168625 


R(l,l) u 


10 


<0.012 




<0.54 




#(1,0) 


10 


<0.012 




<0.33 


A Cep 


#(i,ir 


10 


<0.012 




<0.54 




#(1,0) 


10 


<0.012 




<0.33 



Note. - f lsr is the observed line of sight velocity in the local standard of rest frame. 
FWHM is the line full width at half maximum (for the purpose of calculating column density 
upper limits, the FWHM is assumed to be 10 km s" 1 for all spectra without absorption lines). 
W\ is the equivalent width of the line in Angstroms. a(W\) is the one standard deviation 
uncertainty of the equivalent width. Upper limits for the equivalent width were found by 
taking 3xa(W\). iV(Hj~) is the Bj~ column density. ar(N) is the one standard deviation 
uncertainty of the H^ column density. Upper limits for the column density were found by 
taking 3xa(iV). 



Table 4. Sightline Parameters 





1 v \ ±A S J tot 


a(N) 




E(B - V) 


ri 






L 




Object 


(10 14 cm- 2 ) 


(10 14 cm" 2 ) 


(K) 




(10 21 cm- 2 ) 


f 


(cirr 3 ) 


(vc) 


(10~ 16 


HD 20041 


1.6 


0.18 


76 c 


0.72 d 


4.18 m 


0.67 q 


250* 


5.4 


2.9 


HD 21389 


1.0 


0.08 


51 


0.57 d 


3.31 m 


0.67 q 


250* 


4.3 


1.8 


C Per 


0.7 


0.06 


28 


0.31 d 


1.59 no 


0.60 r 


215 u 


2.4 


3.2 


X Per 


0.8 


0.17 


30 


0.59 d 


2.20 np 


0.76 r 


325 u 


2.2 


3.1 


HD 169454 


0.6 


0.09 


180 c 


1.12 d 


6.50 m 


0.50 s 


265 u 


7.9 


0.9 


HD 229059 


3.9 


0.16 


28 


1.71 d 


9.92 m 


0.67 q 


250* 


13 


2.9 


BD -14 5037 


0.6 


0.16 


26 


1.55° 


8.99 m 


0.67 q 


250* 


12 


0.5 


W40 IRS la 


3.4 


0.37 


27 


2.90 f 


16. 8 m 


0.67 q 


250* 


22 


1.5 


WR 104 


2.3 a 


0.25 a 


38 a 


2.10 s 


12. 2 m 


0.67 q 


250* 


16 


1.4 


WR 118 


6.5 a 


0.18 a 


40 a 


4.13 s 


24.0 m 


0.67 q 


250* 


31 


2.0 


WR 121 


2.2 a 


0.28 a 




1.68 s 


9.74 m 


0.67 q 


250* 


13 


1.7 


Cyg 0B2 12 


3.8 a 


0.36 b 


27 a 


3.35 h 


19.4 m 


0.67 q 


300 v 


21 


1.8 


Cyg 0B2 5 


2.6 a 


0.19 a 


47 a 


1.99 h 


11. 5 m 


0.67 q 


225 v 


17 


1.5 


HD 183143 


2.3 a 


0.08 a 


31 a 


1.28' 


7.42 m 


0.67 q 


250* 


9.6 


2.3 


HD 21483 


<2.2 






0.56 d 


3.25 m 


0.67 q 


250* 


4.2 


<5.7 


40 Per 


<0.9 






0.24J 


1.67 no 


0.35 r 


80 w 


6.7 


<2.6 


o Per 


<0.6 






0.31 d 


1.52 no 


0.54 r 


265 u 


1.9 


<5.0 


e Per 


<0.5 






0.09 j 


0.35 no 


0.19 r 


15 x 


7.5 


<2.4 


i Per 


<0.5 






0.33 d 


1.82 no 


0.38 r 


300 x 


2.0 


<4.5 


62 Tau 


<2.7 






0.37 d 


2.19 n p 


0.56 r 


280 u 


2.5 


<14 


o Sco 


<0.5 






0.73° 


4.23 m 


0.67 q 


225 v 


6.1 


<0.9 


HD 147889 


<0.6 






1.07 d 


6.21 m 


0.67 q 


525 v 


3.8 


<1.6 


COph 


<0.3 






0.32 d 


1.40 no 


0.65 r 


215 u 


2.1 


<1.5 


HD 168625 


<0.8 






1.48 c 


8.58 m 


0.67 q 


250* 


11 


<0.8 


A Cep 


<0.8 






0.57 d 


2.80 np 


0.49 r 


115 u 


7.8 


<1.3 


HD 168607 


<0.6 a 






1.61' 


9.34 m 


0.67 q 


250* 


12 


<0.5 


HD 194279 


<1.2 a 






1.22 1 


7.08 m 


0.67 q 


250* 


9.1 


<1.3 



Table 4 — Continued 







a(N) T(H+) 


E(B - V) 






Tin 


L 




Object 


(10 14 cm- 2 ) 


(10 14 cm- 2 ) (K) 


(mag) 


(10 21 cm- 2 ) 


/ 


(cm -3 ) 


(PC) 


(10- 16 s- 1 ) 




<0.7 a 




0.44 k 


2.55 m 


0.67 q 


250* 


3.3 


<2.1 


P Cyg 


<0.6 a 




0.63 1 


3.65 m 


0.67 q 


250* 


4.7 


<1.2 



Note. — 7V(H+) tot is total Hj" column density. <j(N) is one standard deviation uncertainty of the total 
column density. Upper limits for the column density were found by taking 3xa(N). T(H^) is the excitation 
temperature of fhj~ as determined from the column densities of the (1,0) and (1,1) states. E(B — V) is the 
color excess. is the column density of hydrogen nuclei. / is the molecular hydrogen fraction, rifj is the 
number density of hydrogen nuclei. L is the cloud path length assuming a uniform distribution of gas. £ p is 
the primary cosmic-ray ionization rate. Upper limits on £ p were calculated using 3xer(iV). 

a from lMcCall et all |2002h 



b from lMcCall et all (|l998l ) 

c these high temperatures are most likely caused by inaccurate measurements of the (1,1) state column 
density due to atmospheric interference 

( 2003h 

Thorburn et"all |2003h 



2 derived from method used in 



f derived from lShuping et alJ (jl999n assuming Ry — Ay /E(B — V) = 3.1 
s derived from bendleton et alJ |l994h assuming R v = A V /E(B - V) = 3.1 



h from|Sdmlti dl958h 



CO 
CO 



from lSavage et al.l (|1977l ) 



k derived from intrinsic color of IWegnerl (| 1994T ) 
(|l98i 



from lLamers et al 



•"calculated from N H w E(B - V) x 5.8 x 10 21 cm" 2 mag" 1 in lBohlin et alj (|1978l ) 
"calculated from observed H and H2 column densities 
°from B. Rachford (private communication) 
Pfrom lRachford et all |2002h 

q we adopt / = 0.67 for sightlines without measured column densities 
r H2 fraction derived from same sources as column densities 



s / = 0.5 assumed by ISonnentrucker et al.l (|2007f ) when calculating nn 
* adopted number density 



"from ISonnentrucker et al.l ( 2007 ) 



v from n(H + H2) in ISonnentrucker et al.l (|2007l ) assuming / = 0.67 
w derived from pressure in Ijenkins et all |l983h 
x from|jura| dl975h 



00 
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Table 5. Primary Cosmic-ray Ionization Rate, ( p (10 



-16 S -V 



, for Select Sightlines 



reference 


C Per 


o Per 


e Per 


f Per 


C Oph 


1 


3.2 


<5.0 


<2.4 


<4.5 


<1.5 


2 


0.22 


2.50 


0.01 


0.06 


0.17 


3 


0.17 


1.30 




<0.26 




4 


1-2 


>8 






>4 


5 


5.2 










6 


2.5 











Note. — The upper limits from th is paper are calculat ed using the 3a uncertainty in the 
Hjj" column density. The value from iMcCall et al.l (120031 ) is found by using the conversion 
factor given in equation (flQj) . 
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Federman et al l (119961 ): (4) 



